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The ground and low-lying collective states of a rotating system of A = 3 bosons harmonically 
confined in quasi-two-dimension and interacting via repulsive finite-range Gaussian potential is 
studied in weakly to moderately interacting regime. The A-body Hamiltonian matrix is diagonalized 
in subspaces of quantized total angular momenta 0 < L < 4A to obtain the ground and low-lying 
eigenstates. Our numerical results show that breathing modes with A-body eigenenergy spacing 
of 2hujj_, known to exist in strictly 2D system with zero-range (5-function) interaction potential, 
may as well exist in quasi-2D system with finite-range Gaussian interaction potential. To gain an 
insight into the many-body states, the von Neumann entropy is calculated as a measure of quantum 
correlation and the conditional probability distribution is analyzed for the internal structure of 
the eigenstates. In the rapidly rotating regime the ground state in angular momentum subspaces 
L — IA (A — 1) with 5 = 2, 4 is found to exhibit the anticorrelation structure suggesting that it may 
variationally be described by a Bose-Laughlin like state. We further observe that the first breathing 
mode exhibits features similar to the Bose-Laughlin state in having eigenenergy, von Neumann 
entropy and internal structure independent of interaction for the three-boson system considered 
here. On the contrary, for eigenstates lying between the Bose-Laughlin like ground state and the 
first breathing mode, values of eigenenergy, von Neumann entropy and internal structure are found 
to vary with interaction. 
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I. INTRODUCTION 

Advances made in past few decades have given exper¬ 
imentalists a decisive control over several parameters of 
physical interest of quantum gases such as the density, 
the effective dimensionality, the inter-particle interaction 
0,i and has even made it feasible to confine the desired 
number of particles as in a microchip trap ai- Follow¬ 
ing this, considerable effort has been devoted to explore 
the physics of few-body systems of interacting particles 
a , as it promises to provide a bridge between the micro¬ 
scopic and the macroscopic ensembles a 0 ■ As a con¬ 
sequence, an increasing number of few-body phenomena 
are being studied in a variety of systems under different 
physical situations a^ESl- There are other compelling 
reasons to consider systems with few particles, for in¬ 
stance, a strongly correlated state of the quantum gas is 
experimentally accessible only for small systems, such as 
an analogue of the Laughlin-like state |13l - [l6l| . Another 
reason is that few-body systems allow a higher level of 
control [13, [Hi, to study its role as building blocks of 
strongly correlated quantum many-body states 0 . 

Recent studies have demonstrated that the breath¬ 
ing mode [H [13 is ideally suited to be employed as 
a diagnostic tool for probing the ground and the ex¬ 
cited states of quantum gases in trapped atomic vapors 
[111 [23 paving the way for a novel kind of spectroscopy 
of ultracold trapped gases. For classical system of in- 
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teracting particles confined in 2D by external poten¬ 
tials, breathing mode has been studied earlier [23l - l^ . 
Quantum few-body systems realized in lower dimensions 
with high-precision control over its physical parameters 
serves to motivate the study of breathing mode dynam¬ 
ics [ 1 , [13, llbl - fH - An understanding of few-body sys¬ 
tems may provide an insight into the beyond mean-field 
physics of macroscopic ensembles. 

In this work, we present an exact diagonalization study 
of ground and excited states including the breathing 
modes in a system of three spinless bosons, harmoni¬ 
cally confined and interacting via finite-range Gaussian 
repulsive potential. The results obtained through exact 
diagonalization method attribute a 2fiuj± spacing that 
exists between specific collective excitations, referred to 
as breathing modes of the system in 2D harmonic trap. 
In particular, the paper focuses on the case of rapidly 
rotating bosons (in quantum Hall regime) to explore the 
quantum correlation as well as internal structures of the 
many-body ground and collective excited states. Systems 
in fractional quantum Hall regime characterized by finite 
number of particles or high angular velocity approach¬ 
ing the centrifugal limit of the confining harmonic trap 
are of special interest [H, HI, 113 as the ground state 
of such a system is the well studied Bose-Laughlin state 
[slMUjl - Anticipating experiments with few Bose atoms, 
the purpose of this study is to further extend theoretical 
understanding in this direction. 

The paper is organized as follows. Section|H]provides a 
brief description of the model Hamiltonian for a rotating 
Bose system with hnite-range Gaussian interaction po¬ 
tential and confined in quasi-two-dimensional (quasi-2D) 
harmonic trap. In section imi we present the results for 
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a system of three ultra-cold Bose atoms to investigate 
the physics of breathing modes. For weakly to moder¬ 
ately interacting regime, the low-energy eigenspectra is 
presented to examine the ground and low-lying collec¬ 
tive excitations in rapidly rotating regime. In order to 
gain further insight into the few-body quantum states, 
we obtain the von Neumann entropy as a measure of 
quantum correlation and the conditional probability dis¬ 
tribution to analyze the internal structure (spatial cor¬ 
relation). Section HVl presents summary and outline the 
conclusions of the present work. 


II. THEORETICAL MODEL 
A. The system and the Hamiltonian 

We consider a system of interacting spinless bosons, 
harmonically confined and subjected to an externally im¬ 
pressed rotation with angular velocity fl = Oe^ about the 
z-axis. We assume a stiff confinement of the harmonic 
trap along the axis of rotation so that the axial energy- 
level spacing far exceeds the radial energy-level spacing 
and other energy scales like the interaction energy and 
the rotational energy, yielding an effectively quasi-2D 
system with x-y rotational symmetry. Choosing ?iuj± and 
a_L = ^Jfl/Mu)i_ as units of energy and length respec¬ 
tively, our system of N spinless bosons each with mass 
M and radial confining frequency ujx is described in the 
co-rotating frame by the Hamiltonian — ilL, 

where 


jjlab 



( 1 ) 


Here H = n/uj± (< 1) is the dimensionless angular ve¬ 
locity and L (scaled by h) is the z projection of the to¬ 
tal angular momentum operator. The first two terms in 
the Hamiltonian ([1]) correspond to the kinetic and poten¬ 
tial energies. The third term U (jr^ — rjj) arises from the 
two-body interaction assumed to be Gaussian in particle- 
particle separation [H, [sj 
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with (T_l (scaled by «_!_) being the effective range of two- 
body interaction. The dimensionless parameter g 2 = 
47 ras/a_L measures the strength of interaction where Og 
is the s-wave scattering length for low-energy particle- 
particle collision. 

Recent experimental advancements in atomic physics 
have made it possible to tune Og in ultra-cold atomic 
vapors using Feshbach resonance Qi- Accordingly, in 
a theoretical study, one can vary Og to achieve the de¬ 
sired value of interaction parameter g 2 relevant to the 
trapped model system. We take the scattering length 


to be positive (og > 0) so that the effective interaction 
is repulsive. The finite-range Gaussian interaction po¬ 
tential in Eq. ([2]) is expandable within a finite number of 
single-particle basis functions and hence computationally 
more feasible [3^[3^. In the limit ct_l —?► 0, the normal¬ 
ized Gaussian potential in Eq. smoothly reduces to 
the zero-range contact potential g^J (r,; — rj ) which has 
widely been used in earlier studies [37|. 

The system described by the Hamiltonian in Eq. (HD 
has cylindrical symmetry with respect to the z-axis which 
implies that the z-projection of the total angular mo¬ 
mentum is conserved i.e. L is a good quantum number. 
To obtain the eigenenergies and the corresponding eigen¬ 
states of the N boson system, we employ exact diagonal- 
ization of the Hamiltonian matrix in different subspaces 
of L with inclusion of lower as well as higher Landau lev¬ 
els in constructing the iV-bo^ basis states. 

In Rayleigh-Ritz scheme employed here, the N- 
body variational wavefunction T (ri, r 2 ,..., tat) is con¬ 
structed as linear combination of the symmetrized prod¬ 
ucts (ri, r 2 ,..., Tat)} of a finite number of single¬ 
particle basis functions {Mn.m.n,, (r)}, chosen to be 
the eigenfunctions of the non-interacting single-particle 
Hamiltonian 

^sp = I (-Vi + ri) - H 4 + i (-Vl + Xlz^) , (3) 


identified as the quasi-2D harmonic oscillator Hamil¬ 
tonian in a rotating frame with iz being the single¬ 
particle angular momentum. The eigensolutions of 
HspUn,ni,nA'^) = en.m.n.Mn.m.n, (r), in dimensionless 
form, are known to be: 


= (n -I- 1 — mH) -|- (n^ -|- 1/2) 


M ^ , l {lA-\M})WA/7r^ (,2 

X 4-^ rfl (ri) (A.z^) (4) 


where n = 2nr + \m\ with = 0 , 1 , 2 , • • • and m = 
0, ±1, ±2, • • •. Here lA , ,, (A ) is the associated La- 
guerre polynomial and (Aj;Z^) is the Hermite poly¬ 
nomial. Also Ur = ^{n — \m\) is the radial quantum 
number and m is the single-particle angular momentum 
quantum number. The system here has been assumed to 
be quasi-2D since there is practically no excitation along 
the relatively stiffer z-axis and we, therefore, set Uz = 0 
in Eq. o implying that all the particles occupy only 
the lowest-energy state ^ 0 ( 2 :) = (Az/tt)^/^ of z 

co-ordinate degree of freedom. Therefore Eq. O can be 
written as 


.,m = (n + 1 — mf 2 ) -|- Az /2 , with n = 2nr + 

(r) = . (^rl+X.z-)/2 

i C^{n+\m\})\\ 

X 4-^ rfl (rl) . 


(5) 
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Restricting to = 0 and taking m > 0 in the above 
equation corresponds to the LLL approximation. Taking 
Ur > 0 and allowing m to take positive as well as negative 
values corresponds to going beyond LLLs [H, . The 

iV-body variational wavefunction is 

(ri,r2,... ,rAr) = (ri,r2, • ■ • ,i’Ar) (6) 

where {Ci,} are the variational parameters. The many- 
body index u = (i^oj ■ j • j t'k) labelling the 
many-body basis function (ri, r 2 , ..., rjv) stands for 
a set of single-particle quantum numbers {j = (n,m)} 
and their respective occupancies {iq}. In the present 
work we employ beyond lowest-Landau level approxima¬ 
tion, constructing many-body basis functions {$^,1 us¬ 
ing the single-particle basis Un,m (r) with radial quan¬ 
tum number Ur = ^ {n — \m\) > 0 and angular momen¬ 
tum quantum number \m\ > 0. In the second-quantized 
notation, the Bose field operator can be expanded in 
terms of single-particle basis states as Vj (r) = bjUj (r). 
In occupation-number representation, the A'^-body basis 
function |d>i^) is written in second-quantized form as: 

= n “^(^1) = \vo i^k) (7) 

j=o V J’ 

with X]y=o ^ ^ where j = 

(nj,mj). With these constraints, only the most impor¬ 
tant Fock states from the full basis with a given L (the ac¬ 
tive Fock space) are included. Once the active Fock states 
are constructed as in 0, we diagonalize the Hamiltonian 
matrix. Details of the diagonalization scheme and be¬ 
yond lowest Landau level approximation employed here, 
has been presented in Ref. |^ . 

B. Diagonalization of the Hamiltonian 

Once the active Fock states are constructed, we cal¬ 
culate the matrix elements and subsequently diagonalize 
the Hamiltonian matrix. For N = 5 bosons, we have 
carried out calculations for all the total angular momen¬ 
tum states in the regime 0 < L < AN. Diagonalization 
of the Hamiltonian matrix is performed for each of the 
subspaces of L separately. We have set = 0 in the 
single-particle basis function {z) since there is prac¬ 
tically no excitation along the relatively stiffer z-axis. For 
a given subspace L, the single-particle basis Un,m {r±,(j}) 
spanning the 2D xy plane is chosen as follows. 

It is convenient to define iz = \L/N] where for real 
X the symbol [x] denotes the greatest integer less than 
or equal to x. The single-particle angular momentum 
for the basis functions is now chosen to be: m = — 

rib, (-Z — nb + 1, ■ ■ ■ £z + rib — I, £z + rib, where rib is 
some positive integer that we have chosen to be 3, 4 or 
more depending on the strength of the interaction and 
the computational resources available {rib is a kind of the 


size of the single-particle basis chosen for calculation for a 
given value of T). The single-particle basis functions thus 
chosen are used to construct the variational trial function 
4' = of system for the given value of total 

angular momentum L. 

Since, the system is subjected to an externally im¬ 
pressed rotation about z-axis with angular velocity D, 
we diagonalize the many-body Hamiltonian in given 
subspaces of L to obtain the energy in the corotat¬ 
ing frame E'^°* {L,il,) = [L) — ilL. This can be 

seen as the minimization of E^°‘^{L) subject to the con¬ 
straint that the system has angular momentum expecta¬ 
tion value L and the angular velocity D is then the corre¬ 
sponding Lagrange multiplier. Fixing L fixes D and ac¬ 
cordingly we mention L(f2) instead of rotational angular 
velocity D, in all the tables and figures in the manuscript. 

III. RESULTS AND DISCUSSION 

The results presented here are for a system of A = 3 
Rubidium-87 Bose atoms confined in a quasi-2D har¬ 
monic trap, interacting via repulsive finite-range Gaus¬ 
sian potential. The confining trap frequency is taken to 
be a;_L = 27r x 220 Hz with the z-asymmetry parameter 
\z = uJz/uJx = so that the system has small exten¬ 
sion az = \/fi/Mujz in the z-direction and the dynamics 
along this axis is assumed to be completely frozen. The 
parameters of the two-body interaction potential ([2|) have 
been chosen (t_l = 0.1 and the s-wave scattering length 
in weakly to moderately interacting regime as Og = lOap, 
lOOao, lOOOoo with oq = 0.05292 nm being the Bohr 
atomic radius. The corresponding values of the dimen- 



FIG. 1. The variation of interaction energy Eint contribution 
of the lowest eigenstates with quantized total angular mo¬ 
mentum L {i.e. yrast line), for the system of A = 3 bosons 
interacting via Gaussian potential © with fixed value of inter¬ 
action range ux = 0.1. For a given L, the interaction energy 
is obtained as Eint{L, g2) — E{L, g2 — 0.09151) — E{L, g2 = 
0.0), where E{L, g 2 ) is the total energy of the system includ¬ 
ing the one-body as well as two-body energy terms such as 
kinetic, potential, rotational and interaction energies, in units 
of hijj±. 
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TABLE I. (Color online) The eigenenergy E{Li) of the Li states for N = 3 bosons in the total angular momentum regime 
0 < 1/ < 12 with interaction parameters g 2 = 0.09151 and a± = 0.1 of the Gaussian potential The eigenenergy (in units 
of tioj±) of first ten low-lying eigenstates is shown here. The Li states corresponding to i = 1 (first row in the table) are the 
yrast states (or ground modes) and Li states O2, I2, 23, 84, ds, 65, 67 , 7 ?, 87 , Og, IO9, lls, I29 are the first breathing modes. 
The states O4, le and 28 are the second breathing modes. 


i 

0 

II 

‘Zl 

L = 1 

L = 2 

L = 3 

L = 4 

L = 5 

CD 

II 

>-4 

II 

^-4 

00 

II 

>-4 

L = 9 

0 

II 

^-4 

L = 11 

L = 12 

1 

7.2714 

8.2714 

9.2571 

10.2498 

11.2498 

12.2507 

13.2426 

14.2426 

15.2426 

16.2426 

17.2426 

18.2426 

19.2426 

2 

9.2713 

10.2571 

9.2715 

10.2571 

11.2535 

12.2522 

13.2503 

14.2509 

15.2498 

16.2426 

17.2426 

18.2471 

19.2426 

3 

9.2715 

10.2640 

11.2506 

10.2714 

11.2571 

12.2559 

13.2522 

14.2522 

15.2515 

16.2506 

17.2491 

18.2493 

19.2503 

4 

11.2501 

10.2713 

11.2571 

12.2498 

11.2715 

12.2713 

13.2528 

14.2523 

15.2522 

16.2512 

17.2514 

18.2516 

19.2512 

5 

11.2571 

10.2715 

11.2614 

12.2512 

13.2481 

14.2436 

13.2565 

14.2556 

15.2544 

16.2522 

17.2521 

18.2522 

19.2518 

6 

11.2587 

12.2464 

11.2686 

12.2540 

13.2510 

14.2506 

13.2714 

14.2711 

15.2706 

16.2524 

17.2523 

18.2536 

19.2522 

7 

11.2641 

12.2504 

11.2715 

12.2571 

13.2513 

14.2515 

15.2426 

16.2426 

17.2434 

16.2554 

17.2545 

18.2697 

19.2542 

8 

11.2669 

12.2513 

13.2454 

12.2586 

13.2523 

14.2519 

15.2451 

16.2432 

17.2455 

16.2711 

17.2705 

20.2432 

19.2703 

9 

11.2713 

12.2581 

13.2507 

12.2626 

13.2553 

14.2530 

15.2498 

16.2479 

17.2479 

18.2426 

19.2426 

20.2451 

21.2426 

10 

11.2715 

12.2595 

13.2524 

12.2699 

13.2582 

14.2555 

15.2512 

16.2505 

17.2506 

18.2430 

19.2431 

20.2475 

21.2429 


sionless interaction parameter g 2 = 47 ras/a_L turn out 
to be 0.009151, 0.09151 and 0.9151. The few-body eigen¬ 
states are obtained by diagonalizing the Hamiltonian ma¬ 
trix in each of the subspaces of total angular momentum 
0 < L < AN corresponding to slowly to rapidly rotating 
regime. The low-energy eigenspectra for iV = 3 bosons 
with interaction parameter g 2 = 0.09151 is presented in 
Table U exhibiting how the ground and excited state en¬ 
ergies evolve as L is increased. Eigenstates in a given 
subspace of total angular momentum (columns in Ta¬ 
ble U) constitute a L series (or band). The ith eigenstate 
of the L series is denoted by Li and the corresponding 
eigenenergy by E {Li). The lowest energy eigenstate (cor¬ 
responding to i = 1 ) with angular momentum Li is re¬ 
ferred to as the yrast state of the L series [4l| . The yrast 
line is drawn by plotting the interaction energy contri¬ 
bution of the lowest-energy eigenstate for each of the L 
subspaces [H, In Fig. [H we present interaction en¬ 
ergy (in units of huj±) of the yrast states for total angular 
momenta 0 < L < AN. The red solid circles joined by 
blue line denotes the yrast line. The initial points of the 
plateaus at T = 0, 3 and 6 are the stable ground states 
(other points on the plateaus correspond to metastable 
states). 

Breathing modes. Pitaevskii and Rosch demon¬ 
strated that the Hamiltonian of a harmonically confined 
2D system with zero-range ((5-function) interaction po¬ 
tential possesses SO{2, 1) symmetry due to the scaling 
behavior of (5-function interaction potential. This leads 
to an eigenenergy spectrum with energy spacing of 2huj± 
between two adjacent breathing modes, describing pul¬ 
sation of the system. The above proposition of breathing 
modes for the non-rotating case [ 1 ^ has further been 
generalized to slow rotating as well as rapi(^ ro¬ 

tating regime [i^ . We draw upon these work [1, [23 to 
examine the breathing modes in the following. 

It is observed from Table U that for the non-rotating 
state L = 0 corresponding to H = 0, the energy inter¬ 


val E{02) - E{0i) = 9.2713 - 7.2714 = 1.9999 (in units 
of huj±). The states 0i and O 2 are respectively the yrast 
state and the first breathing mode in the subspace L = 0. 
It is further observed that for the rotating states L > 0 
too, the energy interval has a value close to 2. For exam¬ 
ple, E{ 23 ) - E{2i) = 1.9935 for L = 2, ^( 65 ) - F(5i) = 
1.9929 for L = 5 and ^(lls) - F(lli) = 2.0006 for 
L = 11, all close to 2 (in units of fnjj±). Thus the breath¬ 
ing modes demonstrated to exist in a strictly 2D sys¬ 
tem with zero-range interaction potential are found to be 
observed here in a more realistic quasi-2D system with 
finite-range Gaussian interaction potential ([2]). With an 
aim to study the physics of breathing mode in rapidly 
rotating regime, we focus on T = 6 and L = 12 angular 
momentum states of the rotating three-boson system. 

Breathing modes in rapidly rotating regime. Increase 
in total angular momentum leads to increase in (rota¬ 
tional) kinetic energy and decrease in repulsive interac¬ 
tion energy of the ground state due to centrifugal ac¬ 
tion which moves the bosons apart making their positions 
more correlated. For the ground state of a rapidly rotat¬ 
ing system in 2D with L = N{N — 1), the interaction 
energy reduces to zero and the state is found to be the 
so called strongly correlated Bose-Laughlin state [3l|, 

oc (zi - Zj)‘^ exp (^ 

i<j \ i=l ) 

with L = ^N {N — 1) and q = 2,A,--- ( 8 ) 

Here (zi, Z 2 , • • •, Zat) = {zi} with Zi = ne^^* denote the 
dimensionless co-ordinate of the zth boson in the complex 
plane. For q = 2,4, • • •, the Bose-Laughlin state ([5]) be¬ 
comes an exact eigenstate of the interaction potential in 
Eq. Q with limiting value ct_l —>■ 0 z. e. the 5-function po¬ 
tential. The filling fraction z/, defined (in the thermody¬ 
namic limit) as the ratio of the total number of particles 
to the average number of vortices is given hy u = 1/q for 







5 


the state The exponent q, therefore, 

fixes the filling fraction or equivalently the total angular 
momentum and hence symmetry of the wavefunction. 

In the present study of fV = 3 rapidly rotating bosons, 
we confine ourselves to the first two values q = 2 and 
<7 = 4 for which the filling fractions are i/ = 1/2 and 1/4 
respectively and the corresponding angular momenta are 
L = 6 and L = 12. The yrast states 6 i and 12i in Table H] 
appear as the q = 2 and <7 = 4 Bose-Laughlin states, re¬ 
spectively. 

We first consider the L = 6 subspace to examine the 
ground state and the low-lying excited states, obtained 
variationally through exact diagonalization. In Table |lTl 
we present the eigenenergy E{Li = 6i) as well as the 
von Neumann entropy S{Li = 6 ^) of low-lying Li states 
for three representative values of the repulsive interac¬ 
tion parameter g 2 in the weakly to moderately interact¬ 
ing regime. It is seen from the table that the eigenenergy 
of the yrast state 6 i is independent of the interaction 
parameter and the corresponding wavefunction takes the 
limiting form of Bose-Laughlin state ([ 8 ]) with q = 2 for 
which the interaction energy is zero. It is further ob¬ 
served from the table that the energy of the first breath¬ 
ing mode 67 is also independent of the interaction param¬ 
eter and even as the interaction parameter is varied over 
several orders of magnitude, the ordinal position of the 
first breathing mode 67 with respect to the Bose-Laughlin 
state 61 , in the eigenspectrum, remains unchanged [53|. 


TABLE II. For N — 3 rapidly rotating bosons in total an¬ 
gular momentum subspace L = 6 , values of eigenenergy {E) 
and von Neumann entropy (S) of the ground state and the 
low-lying excited states including the first breathing mode, 
with interaction parameters g 2 = 0.009151, 0.09151, 0.9151 
and range crx = 0.1 of the Gaussian potential The states 
61 and 67 correspond to the ? = 2 Bose-Laughlin state and 
the first breathing mode, respectively. All quantities are di¬ 
mensionless. 


i 

g2 = 0.009151 

gj = 0.09151 

g2 = 0.9151 

E{6i) 

S'(6i) 

E{6i) 

S{6i) 

E{6i) 

S(6i) 

1 

13.2426 

1.5570 

13.2426 

1.5570 

13.2426 

1.5570 

2 

13.2434 

1.6314 

13.2503 

1.6316 

13.3149 

1.6380 

3 

13.2436 

1.2690 

13.2522 

1.2714 

13.3349 

1.3273 

4 

13.2437 

1.5335 

13.2528 

1.5289 

13.3409 

1.5082 

5 

13.2440 

1.6606 

13.2566 

1.6608 

13.3792 

1.6645 

6 

13.2455 

1.5296 

13.2714 

1.5303 

13.5158 

1.5445 

7 

15.2426 

2.3204 

15.2426 

2.3204 

15.2426 

2.3204 

8 

15.2429 

2.2640 

15.2451 

2.2636 

15.2684 

2.2587 

9 

15.2434 

2.1230 

15.2498 

2.1234 

15.3152 

2.1279 

10 

15.2435 

2.2826 

15.2512 

2.2784 

15.3266 

2.2492 


In order to measure the quantum correlation in vari¬ 
ationally obtained states, in particular the breathing 
modes, we calculate the von Neumann (entanglement) 
entropy defined in terms of single-particle reduced den¬ 


sity matrix pi [4^-[5l| 


as 


S'=-Tr(/5ilnpi) (9) 

in subspaces of total angular momentum L. In Table HIl 
we present the von Neumann entropy for N = 3 bosons 
in L = 6 subspace for three different values of interac¬ 
tion parameter. We observe that the value of the von 
Neumann entropy (S) for the first breathing mode 67 is 
large compared to the Bose-Laughlin (ground) state 61 
in L = 6 subspace. Surprisingly, we further note that the 
value of S for both of these states (61 and 67 ) remains 
unchanged as the interaction parameter is varied. It is, 
however, seen from Table El that the eigenenergy and 
the corresponding von Neumann entropy of the eigen¬ 
states lying between the Bose-Laughlin state 61 and the 
first breathing mode 67 change their values as the in¬ 
teraction is varied, in contrast to the ground (z.e. the 
Bose-Laughlin) state and the first breathing modefsj. 

The internal structure (spatial correlation) of a many- 
body state can be analysed by calculating the conditional 
probability distribution (CPD) [H, (H, defined as 


(^1 Ei/j (r - r^) S (ro - r^) 14-) 
(iV-l)E,('I'l^(ro-r,)|4/) 


( 10 ) 


where jd!) is the many-body eigenstate obtained through 
exact diagonalization and ro = (a;o,27o) is the reference 
point (usually chosen to be the position of high density 
for a few-body system like ours). The CPD can be inter¬ 
preted as the probability of a particle being at position r 
under the condition that another one is located at rg. 

In Fig. [21 we present CPD plots of low-lying eigen¬ 
states for = 3 in angular momentum subspace L = 6 
for the interaction parameter g 2 = 0.09151 and cr± = 0.1 

in Eq. 0 with reference point chosen at Tq = ( y/N, 0 ) 


in units of a± in the x-y plane. The ground state 61 
with two equal peaks in CPD, symmetrically placed with 
respect to rg in Fig. 2(a) exhibits a strong anticorrela¬ 
tion (exclusion) structure implying that the probability 
of finding two or more particles at the same position is 
vanishingly small. This underlines the composite fermion 
structure of the Bose-Laughlin state ([ 8 ]). In Fig. |2(b)| for 
the state 62 , the anticorrelation structure vanishes as two 
of the peaks of state 61 , merge around the trap center to 
form a higher peak. The merged peak around the trap 
centre in CPD of state 62 in Fig. |2(b)[ re-distributes it¬ 
self into two distinct unequal peaks in the state 63 in 
Fig. 2(c) where the relatively higher peak shifts away 
from the trap center and a smaller peak appears at the 
opposite end. In the state 64 shown in Fig. |2(d)[ the 
higher peak of state 63 splits into two equal peaks and the 
smaller peak becomes equally prominent with the three 
strongly correlated peaks forming an equilateral triangle. 
We observe that Figs. |2(b)] through 2(e)[ the probability 
of finding one or more particle at the reference position rg 
increases progressively implying increasing tendency to¬ 
wards bosonic correlation. In fact from the CPD plot in 
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(f) 66 (g) 67 BM (h) 68 (i) 69 (j) 610 


FIG. 2. (Color online) Contour plots for conditional probability distribution (CPD) of low-lying eigenstates in angular momen¬ 
tum subspace L = 6 for N = 3 bosons with gj = 0.09151 and = 0.1 in Eq. The plots (isosurface density profiles viewed 
along 2 :-axis) show the probability distribution of finding a particle at position {x, y) when another particle has been fixed at a 
position of relatively high density chosen to be ro = {xo,yo) = (1.732,0), here. Brown-red regions have the highest probability 
density falling off to blue regions of low probability density. 


Fig. |2(f)| for the state 65 with only one peak of about 0.6 
at Tq, it appears as if all the bosons occupy the same po¬ 
sition resulting in a contracted state exhibiting peaked 
bosonic correlation. Interestingly, the state 65 is pre¬ 
cursor to the expanded state 67 identihed as the lowest 
eigenstate of the hrst breathing band and has internal 
structure shown in Fig. |2(g)| with two equal peaks at 
the opposite ends and one right at the trap center. We 
speculate that this first breathing mode 67 possibly has 
structure similar to the state in Eq. (ITT]) . CPD plots in 
Figs. 2(h)|2(j)' which belong to the first breathing band 
Q can interpreted on the same line. 

In order to examine the effect of repulsive interac¬ 
tion on the hrst breathing mode, we present CPD plots 
in Fig. [3] for angular momentum subspace L = Q with 
two values of interaction parameter g 2 = 0.009151 and 
0.9151. First, like the eigenenergy and the von Neumann 
entropy, the internal structure of the ground st ate 61 too 
is in depen dent of interaction as seen in Fig. |3(a)| and 
Fig. |3(b)| This is consistent with the very form of the 
Bose-Laughlin state m which does not allow two par¬ 
ticles being in the same position. We further observe 
that the internal structure of the hrst breathing mode 
67 remains unchanged for the two values of interaction 
parameter in Fig. 3(a) and Fig. |3(b)[ This is corrobo¬ 
rated by invariant values of the von Neumann entropy 
{S = 2.3204), given at top left in each plot. We thus, 
observe that the q = 2 Bose-Laughlin state and the cor¬ 
responding hrst breathing mode in L = 6 subspace are 
similar in having the eigenenergy, the von Neumann en¬ 
tropy and the internal structure independent of interac¬ 
tion. 

The interaction independence of the eigenenergy, the 
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FIG. 3. (Color online) For N = 3 bosons in angular momen¬ 
tum L = 6 subspace, the conditional probability distribution 
(CPD) plots of the g = 2 Bose-Laughlin (ground) state 61 
and the first breathing mode 67 for interaction parameter (a) 
g 2 = 0.009151, (b) gj = 0.9151 with Gaussian width a± — 0.1 
in Eq. ©• In all CPD plots, brown-red regions have the high¬ 
est probability (of finding a particle) falling off to blue regions 
of low probability. 


von Neumann entropy and the internal structure is also 
seen in the g = 4 Bose-Laughlin state and the correspond¬ 
ing hrst breathing mode in angular momentum subspace 
L = 12. The diagonalization is performed for three dif- 




























































7 


TABLE III. For N — 3 rapidly rotating bosons in total an¬ 
gular momentum subspace L = 12, values of eigenenergy (E) 
and von Neumann entropy (S) of the ground state and low- 
lying excited states including the first breathing mode with 
interaction parameters g 2 = 0.009151, 0.09151, 0.9151 and 
range cr± = 0.1 of the Gaussian potential 13). The states 12i 
and 12g correspond to the g = 4 Bose-Laughlin state and the 
first breathing mode, respectively. All quantities are dimen¬ 
sionless. 


i 

g2 = 0.009151 

gj = 0.09151 

g2 = 0.9151 

E(12i) 

S(12i) 

E(12i) 

5(12i) 

E(12i) 

5(12i) 

1 

19.2426 

1.7014 

19.2426 

1.7014 

19.2426 

1.7014 

2 

19.2426 

1.7962 

19.2426 

1.7961 

19.2426 

1.7962 

3 

19.2434 

1.8970 

19.2503 

1.8971 

19.3192 

1.8977 

4 

19.2435 

1.8965 

19.2512 

1.8966 

19.3260 

1.8961 

5 

19.2435 

1.6963 

19.2518 

1.6967 

19.3326 

1.6944 

6 

19.2436 

1.5044 

19.2522 

1.5067 

19.3364 

1.5159 

7 

19.2437 

1.7961 

19.2542 

1.7959 

19.3570 

1.7953 

8 

19.2454 

1.7515 

19.2704 

1.7520 

19.5110 

1.7593 

9 

21.2426 

2.2549 

21.2426 

2.2549 

21.2426 

2.2549 

10 

21.2426 

2.2923 

21.2429 

2.2922 

21.2459 

2.2917 



(a) g 2 = 0.009151 




-3-2-10123 ° 

X 

(b) g2 


FIG. 4. (Color online) For N = 3 bosons with total angular 
momentum L — 12, the conditional probability distribution 
(CPD) of the g = 4 Bose-Laughlin state 12i and the first 
breathing mode 129 for the interaction parameter (a) g 2 = 
0.009151, (b) g 2 = 0.9151 with Gaussian width = 0.1 in 
Eq. ((3)- In all CPD plots, brown-red regions have highest 
probability (of finding a particle) falling off to blue regions 
of low probability. The internal structure as well as the von 
Neumann entropy (marked on the top left corner in each plot) 
of the Bose-Laughlin state 12i and the first breathing mode 
129 are exactly the same even as the interaction parameter is 
varied over two orders of magnitude. 


ferent values of interaction parameter g 2 = 0.009151, 
0.09151 and 0.9151, and the results are presented in Ta- 
ble lllll It is seen from the table that if( 129 ) —if(12i) = 2, 
for the three values of interaction parameter considered. 
The states 12i and 12g are the g = 4 Bose-Laughlin state 
and the hrst breathing mode respectively in L = 12 sub¬ 
space. We further observe that the eigenenergy, the von 
Neumann entropy as well as internal structure (see Fig. 11 
of 12i and 12g states remain unchanged as the interac¬ 
tion parameter is varied over three orders of magnitude. 
However, values of eigenenergy and von Neumann en¬ 
tropy of the eigenstates lying between the Bose-Laughlin 
state 12 i and the first breathing mode I 29 , vary with in¬ 
teraction as is seen from the Table uni 

In order to compare the internal structure of the low- 
lying eigenstates in L = 6 and L = 12 subspaces at a 
given value of interaction parameter g 2 = 0.09151, we 
further analyze CPD plots in Fig. [3 an d Fig. d We 
observe that CPD prohles in Fig. |2(a)| and Fig. 5(a 


of the anticorrelated q = 2 and g = 4 Bose-Laughlin 
states 61 and 12 i, respectively, are similar in having 
two equal peaks equidistant from the reference point 
rg = (1.732,0). It is also observed that the breath¬ 
ing modes 67 and I 29 in Fig. |2(g)] and Fig. 5(i) (where 
the system is assumed to be in the state of expansion) 
have two equal peaks at the opposite ends symmetrically 
placed with respect to the third (higher) peak as well 
as the reference point r p. W e furthe r obse rve that the 
states 06 and 128 in Fig. |2(f)| and Fig. |5(h)[ respectively, 
preceding the breathing modes 67 and I 29 in the respec¬ 
tive eigenspectrum, possess similar internal structures. 
The only peak in their respective CPDs is found coin¬ 
ciding with the reference point rp implying that all the 
bosons are in the same position and the system is in a 
state of contraction. The close analogy between the inter¬ 
nal structures of the low-lying eigenstates in L = 6 and 
L = 12 subspaces suggests that the (radially contracted) 
states 06 and 128 , appear to play an important role in 
the emergence of the (radially expanded) first breathing 
modes 07 and 12g in respective L subspaces. 


Effect of interaction range. Through our exact diag- 
onalization study on N = 3 bosons, we found that the 
first breathing mode in L = 0 subspace is independent 
of the interaction strength g 2 for the interaction range 
cr_L =0.1 in the Gaussian interaction potential We, 
now, fix the value of interaction strength g 2 = 0.09151 
and examine the effect of (t_l on the hrst breathing mode 
in L = 0 subspace as the interaction range is varied over 
0.1 < (Jjl < 0.5. The results obtained for the ground 
and low-lying excited states including the hrst breath¬ 
ing mode are presented in Table IIVI We observe that 
the energy as well as the von Neumann entropy obtained 
for the ground state 0 i and the hrst breathing mode 
remains unchanged for the interaction range a± < 0.3, 
whereas, beyond cr_L = 0.3, both the energy and the cor¬ 
responding von Neumann entropy vary with a±. Thus 
as a± is increased and the Gaussian interaction poten¬ 
tial deviates signihcantly from the i5-function potential. 
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FIG. 5. (Color online) Contour plots for conditional probability distribution (CPD) of various low-lying eigenstates in L = 12 
angular momentum subspace for N = 3 bosons with g 2 = 0.09151 and a± =0.1 in Eq. The reference point ro is located at 
{xo,yo) = (V^, 0) in units of a±. Contour plots depict the isosurface density profiles viewed along the 2 -axis, where brown-red 
regions have the highest probability (of Ending a particle) falling off to blue regions of low probability. 


TABLE IV. For N — 3 rapidly rotating bosons in total angular momentum subspace L = 6 , the values of eigenenergy (E) and 
von Neumann entropy (S) of the ground state and the low-lying excited states including the hrst breathing mode, with fixed 
value of interaction parameter gj = 0.09151 and the interaction range varied over 0.1 < < 0.5 of the Gaussian potential 0 . 

The states 6 i and 67 correspond to the yrast state and the hrst breathing mode, respectively, in L = 6 subspace. 


i 

o-j. 

= 0.1 

cn 

= 0.2 

cu 

= 0.3 

o-j. 

= 0.4 

cu 

= 0.5 

E{6i) 

S{6i) 

E{6i) 

S{6i) 

E{6i) 

S{6i) 

E{6i) 

S(6i) 

E{6i) 

S{6i) 

1 

13.2426 

1.5570 

13.2426 

1.5569 

13.2427 

1.5565 

13.2429 

1.5557 

13.2432 

1.5544 

2 

13.2503 

1.6316 

13.2501 

1.6317 

13.2499 

1.6320 

13.2496 

1.6326 

13.2494 

1.6322 

3 

13.2522 

1.2714 

13.2518 

1.2723 

13.2515 

1.2768 

13.2510 

1.2903 

13.2505 

1.3227 

4 

13.2528 

1.5289 

13.2525 

1.5284 

13.2521 

1.5255 

13.2516 

1.5193 

13.2510 

1.5126 

5 

13.2566 

1.6608 

13.2562 

1.6608 

13.2556 

1.6612 

13.2549 

1.6620 

13.2543 

1.6634 

6 

13.2714 

1.5303 

13.2706 

1.5302 

13.2693 

1.5301 

13.2677 

1.5298 

13.2659 

1.5295 

7 

15.2426 

2.3204 

15.2426 

2.3202 

15.2427 

2.3193 

15.2430 

2.3175 

15.2433 

2.3143 

8 

15.2451 

2.2636 

15.2450 

2.2698 

15.2448 

2.2761 

15.2447 

2.2798 

15.2447 

2.2799 

9 

15.2498 

2.1234 

15.2494 

2.1279 

15.2488 

2.1492 

15.2480 

2.1858 

15.2439 

2.2156 

10 

15.2512 

2.2784 

15.2508 

2.3039 

15.2504 

2.1849 

15.2497 

2.2087 

15.2491 

2.2136 


the energy of the first breathing mode exhibits deviation 
from the value 2hijj±^ for a strictly 2D zero-range inter¬ 
action potential. Such a study has also been presented 
in Ref. for N = 10 bosons in L = 0 non-rotating 
state. Our results, thus, indicate that the first breath¬ 
ing mode remains independent of interaction strength for 
small values of interaction range ct_l (< 0.3) of the Gaus¬ 
sian potential @. 

The independence of the energy of the ground state 
and the first breathing mode for the interaction range 
cr_L < 0.3 of the Gaussian potential (0) in i = 6 and 
L = 12 subspaces for N = 3, can be understood in terms 
of the proposed variational ansatz in Eq. m for the 
ground state and the first breathing mode. However, for 


ct_l > 0.3, the above variational ansatz breaks down. 


IV. SUMMARY AND CONCLUSION 

In conclusion, we have examined the quantum correla¬ 
tion (von Neumann entanglement entropy) and the inter¬ 
nal structure (CPD) of the ground and low-lying excited 
states of a rotating system of three bosons interacting via 
repulsive finite-range Gaussian potential in a quasi-2D 
harmonic trap. The Hamiltonian matrix is diagonalized 
for given subspaces of quantized total angular momenta 
0 < L < AN in weakly to moderately interacting regime 
to obtain the low-lying eigenstates that provide an insight 
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into the evolution of few-body states versus interaction 
and rotation. Our numerical results support the supposi¬ 
tion that breathing modes, known to exist in a purely 2D 
system with zero-range (5-function) interaction potential, 
are also observed in more realistic quasi-2D system with 
finite-range Gaussian interaction potential. 

In the rapidly rotating regime, the evolution of the 
Bose-Laughlin ground state and the low-lying excited 
states with interaction is studied, by fixing the total an¬ 
gular momenta L = ^N{N — 1) = 6 and L = 12 cor¬ 
responding to q = 2 and q = A respectively, for the 
three-boson system. The Bose-Laughlin state which is 
an eigenfunction of the (zero-range) interaction potential 
with eigenvalue zero is strongly correlated and is indeed 
found to have the anticorrelation (exclusion) structure. 
We find that for the three-boson system studied here, 
the Bose-Laughlin state and the first breathing mode 
exhibit similar features such as the interaction indepen¬ 
dence of eigenenergy, von Neumann entropy and inter¬ 
nal structure when the interaction parameter is varied 
over three orders of magnitude. On the other hand, the 
eigenstates lying between the Bose-Laughlin state and 
the first breathing mode, these quantities, namely, the 
eigenenergy, the von Neumann entropy and the internal 


structure are found to vary with interaction. 

Moreover, the eigenstates preceding the first breathing 
mode in L = 6 and L = 12 subspaces have similar inter¬ 
nal structures. The only peak in their respective CPDs 
is found coinciding with the position of the reference par¬ 
ticle (so that the three bosons in the system are at the 
same position). This feature appear to play an important 
role in the emergence of the first breathing mode (or the 
breathing band). This is in contrast to the anticorrelated 
Bose-Laughlin state where the probability of finding two 
or more bosons at the same position is zero. The results 
obtained indicate that the first breathing mode remains 
independent of the interaction strength for small values 
of the interaction range a± (< 0.3) of the Gaussian po¬ 
tential. For cr± > 0.3, the energy as well as the von 
Neumann entropy of both the ground state and the first 
breathing mode vary with a±. Our study demonstrates 
that the von Neumann entropy and the GPD are pow¬ 
erful theoretical tools to gain insight into strongly cor¬ 
related quantum states. We, however, wish to point out 
that some of the results presented here may be the con¬ 
sequence of special geometries like equilateral triangle, 
linear configuration etc., that a three boson system can 
have. 
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teraction independence of values of eigenenergy and von 
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breathing mode 67 is suggestive of writing a variational 
ansatz for these states as 

( 11 ) 

where g = 2, 4, • • •, is an even positive integer (the nor¬ 
malization factor and the exponential factor symmetric 
in {zi} in the above expression have been omitted for the 
ease of writing). The Jastrow prefactor {zi — Zj) 

with angular momentum Ljp = {N — 1) ensures 
interaction independence of the state 4'5 jy(jv-i) 

The Slater-determinant wavefunction for 

composite fermions with angular momentum Lcf = 
{N — 1) is to be constructed from single-particle 
basis states {wn.m}, Eq. ([S]), from the lowest as well as 
higher Landau levels so that ({-^^i}) varia¬ 

tional energy equal to the eigenenergy of the respective 
(the ground or the first breathing mode) states with to¬ 
tal angular momentum L = Ljp -|- Lcf = f Y (Y — 1). 
For the unique choice of single-particle angular momen¬ 
tum states with m = 0,1,..., (Y — 1) in the lowest Lan¬ 
dau levels, the Slater-determinant with q — 2 becomes 
^iw(jv-i) ~ rii<j ~ reduces to the 

Bose-Laughlin ground state wavefunction of Eq. © 
for q — 2. However, for the first breathing mode, a 
prescription to uniquely construct with the 

above constraints on energy, is hard to find and may be 
determined variationally. 


